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An algorithm for obtaining the Taylor coefficients of an expansion of Feynman diagrams 
is proposed. It is based on recurrence relations which can be applied to the propagator 
as well as to the vertex diagrams. As an application, several coefficients of the Taylor 
series expansion for the two-loop non-planar vertex and two-loop propagator diagrams 
are calculated. The results of the numerical evaluation of these diagrams using conformal 
mapping and Pade approximants are given. 



1. Introduction 

Recently a new method to calculate Eeynman diagrams was proposed^. It is 
based on the small momentum expansionn, conformal mapping and construction 
of Pade approximants from several terms in the Taylor series of the diagram. The 
method was successfully applied to the evaluation of two- and three- loop diagramsEl. 
As it was observed, a suitably accurate approximations to the integrals can be 
obtained with 20-30 coefficients in the Taylor series. The computation of two-loop 
vertex diagrams reveals the necessity of an efficient algorithm for the expansion of 
the diagrams w.r.t. external momenta. 

To outline the problem, let us shortly describe the existing approach to the 
small momentum expansion. At present the only method of expansion is based on 
differentiation of the diagram w.r.t. external momenta. For the propagator type 
integrals the prescription for the small momentum expansion was formulated in 
Rcffl. Any coefficient in the Taylor series w.r.t. the external momentum q 2 can 
be obtained by applying □ = g^-jj^jr in an appropriate power to the diagram and 
putting the external momentum to zero. 
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The general expansion of a scalar 3-point function C(qi,q2) can be written as 



C( gi ,q 2 )= ]T a lmn {ql)\ q 2 2 r(q iq2 r = E E *lmn{ql)\ql) m {<liq2) n . 

/,m,n— L— Z+m+n— L 

(i) 

The coefficients <2; mn are to be determined from a given diagram. They can be 
obtained by applying the differential operators Dy = j^— g^p- several times to both 
sides of ([I]). This procedure results in a system of linear equations for a\ mn . At a 
fixed L, one obtains several systems of [L/2] + 1 equations. 

Differential operators (Df's), projecting out the coefficients aoon in an arbitrary 
space-time dimension d, are: 

Dfo0n - ^ 2T(i)T(n - 2z + 3)I> + d 2)T(n + d/2) ( ° 12) (DllD22) ' 

(2) 

Applying -D/oon to C(qi, 92) and putting the external momenta equal to zero yields 
the expansion coefficients aoo„. For the coefficients a/on the following projection 
operator was obtained: 

r(f + n) 



Dfl0n - r(i + i)r(l + i + n) V—J Dfoan - (3) 

The projection operator for arbitrary ai mn is yet unknown. 

An essential element of the above expansion is multiple differentiation w.r.t. 
multidimensional vectors, ft turns out that the computer implementation of algo- 
rithms for the analytic calculation of the Taylor coefficients based on differentiation 
is not very effective. The computation of multiple sums with multiple differentia- 
tions in an arbitrary space time dimension d, even with the advanced possibilities 
existing in FORM leads to significant computational difficulties. 

2. Recurrence relations 

In the present paper we propose a new approach to this problem, which we 
expect to be more suitable for the computer implementation. The algorithm for the 
small momentum expansion of scalar integrals can be formulated as follows: 

• Firstly, the propagators with external momenta are to be expanded as 

1 _ 1 ga+j-i)! (2k iqi y 



(k-(-2k iqi + q (- m {y (fc?+«?-mf)i^ C?-l)ffl (kf + qf-miy 

This expansion is valid in the whole integration region since after the Wick 
rotation the expansion parameter is always < f . 



3 



• Secondly, by using recurrence relations the external momenta are to be fac- 
tored out. 

• Thirdly, the remaining factors should be expanded w.r.t. the external mo- 
menta qf. 

• Finally, using another kind of recurrence relations, bubble integrals must be 
reduced to a set of master integrals and trivial ones. 

The proposed algorithm does not need any explicit differentiation or solution of 
linear systems of equations. All required operations are efficiently implemented in 
many computer algebra systems. A key element of the algorithm is the recurrence 
relations allowing one to factorize external and integration momenta. In the present 
paper such recurrence relations will be presented for the two-loop case. It's not 
difficult to get them for higher loop diagrams. For the three-loop case they will be 
presented in a future publication. 

We shall now describe a realization of the algorithm for the two-loop propagator 
and 3-point vertex diagrams. At the one-loop level the small momentum expansion 
is more simple, and the algorithm will be evident from the two-loop consideration. 
At the two- loop level there are two topologically different 3-point vertex diagrams 
and one propagator type topology (see Fig. I). We assume that the lines in Fig.l 
correspond to arbitrary powers of the propagators. 
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Both vertex diagrams, after expansion of propagators w.r.t. scalar products 
(hqj), yield integrals of the type: 



d d k! d d k 2 f(ki,k2,qi,q 2 2 )(kiqiy i (k 1 q2y 2 (k2qiy 3 (k2q2y 4 =v(j 1 ,j2,h,ji). (4) 
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These satisfy the following recurrence relation: 



(d + ji + h + h + k - 2 ) (d + ji + h - 3)w(ji , 32 , 33 ,34) = 

{{d + h+h+k ~ 3) [(ii - 1)*? q\ 1~ +hk{ ( qi q 2 ) 2 + j 3 (feifea) q\ 3-] 

- j33i qjki2+3-4- - j 4 (j 4 - 1)^1(91 g2)2 + 4-4- - j 4 (ji - l)(fcifc 2 )<Zi 12+4 

+k{d + h -32+33+34, - 4) (fcifc 2 ) (qiq 2 ) 4"} 1~ <g> j 2 , ja, k), (5) 

where 1 v(ji, ...) = w(ji±l, ...) etc., and the sign <S> means that the scalar products 
of the integration momenta in braces must be considered under the integral sign in 
v. Relation (||) can be derived from the following tensor formula 



after contraction with q±,q2- Here C^{x) are Gegenbauer polynomials, fcifc 2 = 
{k\ki) I \/ k\k\, S is the totally symmetric sum of products of g's and P v means sym- 
metrization w.r.t. the indices v, i.e. P u g v ^S [ -^ ] = {g v ^ S [ -' V3] + g v ^S [ -' V2] + 
g V2U3 S^^' vi]> ) /i, etc. The integral is equal to zero if L + N is odd. In the above 
formula L < N is assumed without loss of generality. 

Expansion of the self-energy diagrams will give integrals which correspond to 
(|J) with 31 = 33 = 0. In this case the recurrence relation is very simple: 



{d+j 2 +k-2)v(0,j2,0,k) --Z 2 2 {02-l)fci2-+j4(fcifc2)4-}2-® V (0,j 2 ,0,j 4 ). (7) 



In RefJ an explicit formula for v(0, 32, 0,^/4) in terms of a onefold sum is given. 
The computer implementation of that formula, however, is less effective than the 
application of the recurrence relation (^) . Substitution of the sum blows up expres- 
sions, producing too many similar terms at once. 

A recurrence relation for the one-loop vertex diagrams is also a special case of 



(§) with 33 = ji = 0: 

(d+ji+j2-2)«(?i,i2,0,0) = fc?[(3i-l)g?l-+j 2 ( qi q 2 ) 2-}l-®v{ji,j2,0,0). (8) 
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Here it is assumed that the integration w.r.t. &2 in v is omitted. 

The recurrence relation (||) can also be used for evaluating anomalous dimensions 
of the moments of structure functions in deep inelastic scattering in the MOM renor- 
malization scheme. In this case, one needs to calculate propagator type integrals 
with vertex insertions containing some additional vector £. The small momentum 
expansion of these diagrams will produce integrals of the type 



d d kxd d kif{kxM,(ll){kxqi) h {kiO h (Mi)'' 3 , (9) 
again corresponding to the special case of (JsJ) with j'4 = 0: 

(d + 31 + 32 + 33 - 2)»(j'i. h, 33, 0) = (10) 
[(ii - l)klq\l- +32^(^1)2- + j 3 (fci*2)g?3-]l- <g> v{j u j 2 , j 3 ,0). 

After some modifications, the recurrence relation (|^) can be used for the calcu- 
lation of the Taylor coefficients in the expansion of diagrams in axial gauges. 

By using relation (|s|) and relations derived from it by exchanging q\ <-» q 2 and 
k\ «-> k2, the exponents of the scalar products can be reduced to zero. One should 
choose the scalar product with the smallest exponent, apply the recurrence relation 
until this exponent is reduced to zero and repeat this procedure until only one scalar 
product in some power remains. For these integrals the following explicit formula 
can then be used: 

f{k x ,k 2 ,qi,ql) {k 1 q 1 f^d d k l d d k 2 (11) 



|V J f(kt, k2,qlql)(kl)^d d k 1 d d k 2 . 



(d/2) n 

After that we obtain bubble-like integrals depending only on the external mo- 
menta squared: 

' d% d d k 2 /(fci,fe2,4i 2 ,<zf) (fc?r(fcife) J2 (fc 2 2 ) J3 . (12) 



The most optimal way to evaluate these integrals depends on the mass values. If 
some masses are the same, or one mass is zero, then it is worthwhile to expand the 
integrand w.r.t. q\,q\ and after that to calculate the bubble integrals using recur- 
rence relations given in Refsflfl. If the masses are different it will be more efficient 
to evaluate the integrals ( |T^ ) using the recurrence relations before the expansion 
in q\^q\ is performed. In this case the presence of q\,q\ in the denominators can 
be considered as a mass shift, and therefore, the same recurrence relations can be 
applied. 

3. Applications 

To check our algorithm, we repeated the calculation of the diagram that occurs 
in the process H — > 77 u. In comparison with the old algorithm based on differ- 
entiation, the execution time reduced by a factor of about 20. For example, with 
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the old algorithm the evaluation of the 28-th coefficient in the Taylor series for this 
diagram took5 hours CPU time on the DEC 3000. With the new algorithm, run- 
ning FORM □ on the same computer, it took 15 min to obtain the result. Further 
improvement in the efficiency of the computation can be achieved if the algorithm 
can be implemented in a multiple precision FORTRAN program 0. In this case, as 
it was observed in Ua , the numerical approach should substantially more efficient 
than the analytical one. 

The method proposed here was used for the evaluation of the Taylor coefficients 
of the nonplanar scalar vertex and propagator diagrams. In both cases all inter- 
nal lines were considered to be massive. For simplicity the same mass was taken 
and q\ — q 2 — was chosen, so that the diagram depended on just one variable 
q 2 = — 2gi(j2. The evaluation of the non-planar diagram by numerical methods was 
complicated due to the Lntegrable singularities i. It is interesting to investigate how 
the method proposed iiB works in this case. The diagram has a cut for q 2 > 4m 2 . 
After the conformal mapping we find a good convergence of the sequence of Pade 
approximants derived from several first Taylor coefficients. Table 1 demonstrates 
the quality of the convergence and contains values (up to a factor 10 _9 /(167r 2 ) 2 ) 
for the diagram on the cut. Errors were estimated by comparison of the [8/8] and 
[9/9] Pade approximants. With 18 coefficients, as one can see from Table 1, the 



Table 1: Results for [9/9] on the cut {q 2 > 4m 2 ), to = m t = 150 GeV 
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accuracy of calculations for q 2 < 400m 2 is better than 1%. The convergence of the 
approximants below the cut is substantially better than on the cut. In Table 1 we 
also give results on the cut for the propagator diagram. The precision of the results 
for the propagator diagram is practically the same as for the vertex one. This fact 
can be an indication that for the diagrams with arbitrary scalar products of the 
external momenta, taking in (1) L = 18 — 20, we can achieve a similar accuracy. 
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